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ABSTRACT 

Context. Most of our knowledge about star formation is based on studies of low-mass stars, whereas very little is known about the 
properties of the circumstellar material around young and embedded intermediate-mass T Tauri stars (IMTTSs) mostly because they 
are rare, typically more distant than their lower mass counterparts, and their nearby circumstellar surroundings are usually hidden 
from us. 

Aims. We present an analysis of the excitation and accretion properties of the young IMTTS DK Cha. The nearly face-on configuration 
of this source allows us to have direct access to the star-disk system through the excavated envelope and outflow cavity. 
Methods. Based on low-resolution optical and infrared spectroscopy obtained with SofI and EFOSC2 on the NTT we derive the 
spectrum of DK Cha from ~0.6ytim to ~2.5yum. From the detected lines we probe the conditions of the gas that emits the Hi IR 
emission lines and obtain insights into the origin of the other permitted emission lines. In addition, we derive the mass accretion rate 
{Mace) from the relationships that connect the luminosity of the Bry and Pa/3 lines with the accretion luminosity (Lore). 
Results. The observed optical/IR spectrum is extremely rich in forbidden and permitted atomic and molecular emission lines, which 
makes this source similar to very active low-mass T Tauri stars. Some of the permitted emission lines are identified as being excited by 
fluorescence. We derive Brackett decrements and compare them with different excitation mechanisms. The Pa/^/Bry ratio is consistent 
with optically thick emission in LTE at a temperature of ~ 3500 K, originated from a compact region of ~ 5 Rq in size: but the line 
opacity decreases in the Br lines for high quantum numbers n,,,,. A good fit to the data is obtained assuming an expanding gas in LTE, 
with an electron density at the wind base of ~10" cm"^. In addition, we find that the observed Brackett ratios are very similar to those 
reported in previous studies of low-mass CTTSs and Class I sources, indicating that these ratios are not dependent on masses and 
ages. Finally, Lore ~9Lo and Muc-c ~3xlO"' Mq yr"' values were found. When comparing the derived M„cc value with that found in 
Class I and IMTTSs of roughly the same mass, we found that M^cc in DK Cha is lower than that found in Class I sources but higher 
than that found in IMTTSs. This agrees with DK Cha being in an evolutionary transition phase between a Class I and II source. 

Key words, stars: formation - stars:circumstellar matter - stars: pre-main sequence - ISM:individual objects: DK Cha, IRAS 12496- 
7650 - Infrared: ISM 

1. Introduction ing the emission from each of the components of this system 

_ . , ^ can only be probed through spectroscopic studies that allow us 

Intermediate-mass TTaun stars are the young precursors of ^^^^^ quantitative information on processes that take place 

Herbig AeBe stars with special type s rangmg froni K to late .^^^ unresolved regions (such as the accretion through 

F, and ranging from 1 Mo to 5 (e.g.| Ca vet et al|20C)4 . Their ^^^^^^^j ^^^^ ^.^^^.^^ ^^^^^^ ^^^^ .^^^^ ^.^j^^ 

disks dissipate m ~5 Myr, i.e., more quickly than their low-mass .^^^^^^^^ .^^ properties can be studied through forbid- 

H classical TTaun (CTT) counterparts, which makes it very dif- ^^.^^.^^ ^.^^^^ ^^^^ ^^^^.^^^^^^ properties are usually 

, ficult to study them in the early evolutionary phase. Besides, ^^^^^ ^j^^.^^ transitions. In this context, infrared spec- 

the complex structure of the emission region surrounding pro- ^ ^^^.^^ ^ ^^^^^j ^^^j ^^^^ embedded 

tostars makes It more difficult to probe the disk accretion prop- ysOs where the high extinction usually prevents us from ob- 

erties, a crucial step towards understanding disk evolution in ^^^^.^^ ^^^^ ^^^^^^^ ^^^^^^^ ^p^^^^j ^^^^^^^ However, in a 

intermediate- niass stars. , . , , small number of cases, when the object is observed face-on, we 

In the nearby surroundings of protostars one does indeed ^^^^ ^^^^^ ^^^^^^ ^ ^^^^^^ ^^^^^^ 

expect emission from the heated dense envelope, the disk the excavated envelope even at shorter wavelengths. 

jet and the UV-heated and shocked cavity walls, excavated by ^j^^ ^^^^3 j^^^ (M.=2Mo, U„=29.4Lo) is one of 

the outflow. Information about the dmerent phenomena involv- ^, i u »u i » if c 

^ these rare examples where, thanks to its nearly face-on conngu- 

Send offprint requests to: R. Garcia Lopez, e-mail: 
rgarcia@mpifr.de 

* Based on observations collected at the European Southern 
Observatory, La Silla and Paranal, Chile (ESO program 082.C-0264(A) young object still surrounded by its envelope, as shown by the 
and 084.C-0308(A)). large amount of warm molecular gas detected on-source and 
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Fig. 1. Normalised to the continuum EFOSC2+SofI spectra from 0.6 to 2.4 jum of the PMS star DKCha. The strongest Hnes are 
identified. 



analysis and discuss our observations. Finally, we present the 
conclusions in Sect.|5] 



2. Observations and data reduction 

As part of the POISSON project, we observed DKCha in the 
wavelength range from ~0.6 fim to ~2.5 fim using the optical and 
infrare d spectrographs EFOSC2 and SofI at t he ESO-NTT tele- 
scope (iBuzzoni et al.lll984l : lMoorwoo d et sT. 1998). The similar 
spectral resolution and the selection of comparable slit widths 
for both instruments (R~700-900 and shts widths of OfH and 
0'.'6 for EFOSC2 and SOFI, respectively) allowed us to acquire 
optical and infrared spectra of homogeneous characteristics. 

The EFOSC2 spectrum (hereafter "optical spectrum") was 
acquired with grisml6 and roughly ranges from 0.6fj.m to 
1.0 jum. The infrared SOFI spectra were taken using both the 
blue- (GB: 0.95yum-1.64jum) and red-grism (GR: 1.53yum- 
2.52yum). The optical and near-IR (NIR) data were acquired 
close in time on 12''' and 16''' February 2009, respectively. The 
total integration times on DK Cha were 300 s for the optical and 
60 s and 20 s for the blue and red IR grim spectra. These were 
acquired through standard A-B nodding. 

For data reduction standard IRAlQ tasks were used. 
Wavelength calibrations were performed using argon-xenon 

' IRAF (Image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are operated 



the presenc e of a strong molecular outflow (van Kemp en et al] 
l2009l l2006l). In addition, its spectral energy distribution (SED) 
corresponds to that of an YSO in a transition phase from a 
Class I to a Class II source (I van Kempen et al.l2010l : ISpezzi et af] 
[2OO8). These characteristics together with its nearby location 
(d~178pc) make DK Cha an ideal candidate to probe the first 
stages of the formation of Herbig Ae stars. 

Most of our knowledge about this source is based on far-IR 
and sub-millimeter observ ations (e.g.. Ivan Kempen etal1l2006l 
I20 1 Ot iGiannini et aPl 1 9991) . Optical spectr oscop ic studies of DK 
Cha have been made by [Hughes et"al] (1199 li) . In that work 
the first extensive study of the optical and infrared photomet- 
ric variability of DK Cha was also reported, pointing out the 
highly variable nature of this source. Here, we will present 
an optical/IR spectroscopic study of DK Cha aimed at obtain- 
ing insights into the characteristics of its circumstellar material 
and accretion properties. This study was conducted as part of 
the POISSON (Protostellar Objects IR-optical Spectral Survey 
On NTT) project, whose results will be presented in a series 
of papers on the properties of YSOs b elonging to different 
clouds (e.g. Chal-II regions and L1641 lAntoniucci et al.ll201 ll 
ICaratti o (3aratti et al.ll20Ilh . 

In Sect.|2|the observations and the data reduction will be de- 
scribed, while in Sect. [3] the first complete spectrum of DKCha 
from 0.6 to 2.4/im will be presented. In Sect.|4| we show the 
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Table 1. Detected emission lines in DKCha. Table 1. Continued. 
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Notes, (m) indicates lines belonging to multiplets. Observed fluxes, not Notes, (m) indicates lines belonging to multiplets. Observed fluxes, not 
corrected for extinction. corrected for extinction. 



lamps. The NIR spectra were corrected for the atmospheric spec- 
tral response by dividing them by the spectrum of a suitable stan- 
dard star. The three overlapping segments were then joined to 
obtain a single optical/IR spectrum of DK Cha. Unfortunately, 
during the two observing nights variable seeing conditions (see- 
ing in the range of 0.75"-1.25" and 1.0"-2.5" for the optical and 
infrared data, respectively) prevented us from performing pre- 
cise standard flux calibrations using spectro-photometric stan- 
dards. This is because of considerable flux loss owing to the fi- 
nite slit widths. To flux-calibrate the optical and infrared spec- 
tra we therefore employed the acquisition image (taken with 
the narrow-band filter NB 2.191) using four field stars as "stan- 
dard stars". Moreover, additional photometric observations of 
DKCha were acquired eight months later using the infrared 
spectrograph and camera ISAAC on the ESO-VLT telescope 
to test the goodness of the acquisition image photometry. The 
ISAAC naiTow-band filter at 1 .2 1 ^va (where no significant spec- 
tral features were expected to contribute to the total flux) was 
chosen instead of the standard J broadband filter because of the 
bright J-magnitude of DKCha (2MASS J~9.3 mag) that would 
otherwise saturate the detector The photometric data and the ac- 
quisition image were reduced using standard IRAF ' routines. 



by AURA, Inc., cooperative agreement with the National Science 
Foundation. 



The magnitudes derived from both the ISAAC and acquisition 
image photometry agree within 0. 1 mag, indicating no signifi- 
cant variations in the time elapsed between the observations. 

3. Results 

3.1. Detected lines 

The optical and NIR spectra of DK Cha normalised to the stel- 
lar continuum are shown in Fig.[T] The spectra show numer- 
ous emission lines, most of them detected at S/N greater than 
5, whereas no significant absorption features are detected. Line 
identification together with fluxes and equivalent widths are 
given in Table[T] Identification of these lines was made using 
the Atomic Line List database, while the goodness of the classi- 
fication was checked through the detection of several transitions 
from the same atom and multiplet. In addition, the spectra of 
other YSOs ( Kellv et al. 1994) and obje cts w ith expected similar 
excitation conditions (IWalmslev et al.l2000l) were also employed 
to recognise the features. 

The DK Cha spectra show permitted emission lines simi- 
lar to those observed in active YSOs, as for example the pres- 
ence of Ha, He i, the Ca II infrared triplet, and the Pa /3 and Br y 
emission lines. Indeed, the most abundant features observed in 
our spectra are H i emission lines coming from the Brackett se- 
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Table 2. DK Cha photometry. 
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rie. Permitted emission lines from species such as Nai, Mgi, 
C I, O I, and N i are also detected, showing upper state excita- 
tions above SOOOOcm"', 47000cm-', 69000cm-', SSOOOcm"' 
and lOOOOOcm respectively. Several of these lines are usu- 
ally observed in the optical and infrared spectra of CTTSs and 
low-mass Class I sources, as well as Herbig AeBe stars, while 
other lines (such as the Ni lines) may be excited by fluores- 
cence, which we will discuss in Sect. 14. 31 and so metimes associ- 
ated with photon dominated regions (PDRs; e.g. lMarconi et alj 
ll998VWalmslev e t al.ll2000h . 

Forbidden emission lines such as [Feii], [S ii] and [Oi] tran- 
sitions are also detected in our spectra. These lines are com- 
monly accepted to be associated with jets from young stars and 
are often observed in the spectra of low-ma ss CTTSs and Class I 
sources s howing outflow activity (see, e.g. lHartigan et alJll995t 
iNisini et a l. 2005b; Podio et al. 2006). 

Finally, in addition to the atomic lines, emission from mo- 
lecular transitions are also observed in the infrared spectrum of 
DK Cha. In particular, weak H2 emission at 2.12 /urn and band- 
head CO emission around 2.3 fim have been detected. 

3.2. Photometry 

Table|2] shows the photometric value retrieved from the ISAAC 
photometry at 1.21 yum together with the V, and K magni- 
tudes derived from the measured fluxes in the DK Cha spec- 
trum at each eff'ective wavelength. A s already mentioned , DK 
Cha is known to be a variable star. Hughes et al.l (Il991b and 
iMolinari et al.l (Il993h studied the variability of DK Cha in a five- 
year period (from 1987 to 1992). They reported a AJ variability 
up to ~3 mag in 5 years. The K- and J-magnitudes re ported in the 
literature vary between 5- 6 mag and 9-11 mag ('e.g.. lAlcala et al.l 
I200 8'. Hughe s et al.lll99ll) . respectively, w hile the V-magnitude 
ranges from V~ 18.7 m ag dHughes et alj[l991) to V~16.7mag 
(IZacharias et al.l2004l) . Our measured magnitudes are within this 
range of values, indicating that at the time of the observations 
DK Cha was close to its brightness minimum. 

4. Analysis and discussion 

4.1. Extinction 

Deriving the correct extinction value is crucial for the com- 
putation of YSO physical parameters. Frequently, extinction 
values are derived from colour excesses once a spectral type 
and an extinction law are assumed. The visual extinction (Ay) 
derived in this way depends on the adopted extinction law. 
Here, we have derived Ay fro m (I-R), (I-J) and (I - H) as- 
suming the stan dard colours of Kenvon & Hartmamil (Il995h 
fo r a FOV star (ISpezzi et al. I I2008h and the reddening law 
of ICardelU et"aD (!l989^ with Ry=5.5 because the interstellar 
medium in Chameleon clouds shows Ry v alues around 5-6 in 
the densest parts of the cloud (see, e.g. ICovino et al. I CM 
ISpezzi et al.l l2008l) . The discrepancies between the Ay val- 



ues derived from different colours are not larger than ~2 mag 
(Ay=10.4mag, 10.5 mag and 12.6mag for the (I-R), (I-J) and 
(I-H) colours), with an average value of 11 + 1 mag. For these 
estimates no veiling correction to the colours was considered. 
Nevertheless, ignoring the ve iling will not lead to an Ay varia- 
tion greater than 0.5 mag (see lCieza et al.ll2005l) . which is within 
our Ay uncertainty of +1 mag. In addition, the derived Ay value 
agrees well with the extinction derived from the silicate ab- 
sorptio n feature at 10 nm as d e tected in the ISO-SWS spec- 
tra by lAcke & van den Anckej (l2004h . From their spectra of 
DK Cha, an Ay value ranging from ~10mag to ~12m ag was 



derived, depending on whethe r the Riek e & Lebofskvl (I1985L 
Ay/r(9.7)=16.6) or flielM^H (11998, Ay/T(9.7)=19.3) relation- 
ships were used. In the following we adopt an average Ay value 
of 1 1 mag. 

4.2. Hi line emission 

The line ratios among Brackett and Paschen lines can be used 
to obtain information on the physical conditions of the emitting 
gas, which helps us to determine the origin of the Hi lines. The 
SOFI spectra taken with the blue- and red-grism allow us to si- 
multaneously detect a series of Brackett lines from Br y to Br 20, 
and the Paschen series lines Pa/?, Pay and PaS. This allows us to 
minimise errors in the line ratios caused by line flux variabil- 
ity. In addition, line ratios are independent of uncertainties in the 
absolute flux calibration. 

In Fig.|2]the ratios of the different Brackett and Paschen lines 
with respect to the Br -y and Pa /? line (left and right panel, re- 
spectively) are shown as a function of their upper quantum num- 
ber (n,,^). Only unblended lines and those in clear atmospheric 
windows were taken into account. The fluxes of the lines were 
derived from the observed equivalent widths (EWs). Because 
DK Cha is of spectral type FO (Spezzi et al. 2008), the contribu- 
tion of the intrinsic photospheric component to the observed H i 
EWs cannot be neglected. Therefore, the EWs were corrected 
from the intrinsic photospheri c absorption contribution f ollow- 
ing the same procedure as in ' Garcia Lopez et af] (|2006). This 
procedure assumes that the observed flux is the sum of the emis- 
sion from the stellar photosphere, the circumstellar gas and the 
disk. In addition, it assumes that the star dominates the V-band 
emission (see, Rodgers 2001). The latter assumption is further- 
more supported by the fact that IMTTSs are not veiled around 
5500 A (Calvet et al.ll2004.) . The EW of the photospheric com- 
ponent of each H i Une, corrected for the H- and K-band veil- 
ing (r^ ~0 .6, Th ~0.3), was computed using a FO spectral tem- 
plate from lRavner et al . (20090 No veiling was detected in the 
J-band. All relevant parameters used to derive the photospheric 
equivalent widths as well as the stellar parameters are given in 
Tables|2] and |4] The equivalent widths of the circumstellar H i 
line components and the computed fluxes are shown in Table[3] 
The errors bars shown in Fig.|2]coiTespond to an uncertainty of 
+ 1 A in the EWcc. In the same figure, the decrement observed 
in other YSOs is also plotted for comparison: namely, the line 
ratios as observed in the low-mass Cl ass I source HHIOO-IRS 
(IPodio et al ]|2008; Nisini et al. 2 005ah and the average values 
inferred from ^Barv et al..(2008) in a sample of CTTSs located in 
the Taurus cloud. 

The DK Cha Brackett and Paschen ratios are very similar to 
those of low-mass CTTSs and Class I sources, although the T 



' Online at the IRTF spectral library Web 
http://irtfweb.ifa.ha waii.edu/~spex/IRTF_Spectral-Library/| 
Data_Format.html 



site: 
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Fig. 2. Left panel: Ratio of the DK Cha Brackett lines with respect to Br y plotted as a function o f the upper quantu m number 
(n„p; filled black circles). For comparison the Brackett decrements for the Class I source HHIOO-IRS (Nisini etalJliOOl open blue 
triangles) and the weighted mean line ratios for the Brackett series on the CTTSs sample of Barv et al. (2008, open green squares) 
are also included in the plot. Right panel: Ratio of the DK Cha Paschen lines with respect to the Pa/3 line as a function of the 
upper qua ntum num ber n„„ (black filled circles). For comparison, the weighted mean line ratios for the Paschen series on the CTTSs 
sample of lBarv eF al. (2008, open green squares) are also shown. 




Fig. 3. Left panel: Comparison between the observed DK Cha Brackett decrement and the expected curves for Case B recombi- 
nation and black-body emission. The crosses represent the ratios from case B recombination at T=10000K, n=10'"cm"^ (black 
continuous line) and T=1000K, n=10"^ cm"^ (green dashed line), while dots represent the black-body curves at T=3500K (dashed 
red line). Right pa nel: Com parison between the observed DK Cha Paschen decrement (black dots), the Paschen decrement of the 
CTTSs sample of iBarv et alJ (l2008l open green squares) and the expected Paschen ratios from Case B recombination at T=1000K, 
n=10'"cm"^ (green dashed line). 



Tauri sample of iBarv et alJ (l2008l) have slightly lower Brackett 
ratios at N„;, > 13 with respect to DK Cha and HHIOOIRS. This 
might indicate that the excitation conditions giving rise to the 
Brackett series have a weak dependence on the mass and evolu- 
tionary stage of the source. 

4.2.1. Case B recombination 

To constrain the physical conditions where the Brackett and 
Paschen lines originate, we have explo red different excitation 
mechanisms. In particular, iBarv et al] (12008) found that the 
Brackett and Paschen decrements of their sample of CTTSs 



were always fairly well fitted by a standard Case B re- 
combination at low gas temperature (<2000K). In Fig. [3] we 
plot the expected ratios under Case B conditions that fit the 
Bary's sample (n=10'°cm ~^, and T=1000K), derived from the 
iHummer & StorevI (Il987h calculations^. We find that this model 
is unable to reproduce the observed Brackett decrement in DK 
Cha, bacause the line ratios at N,,^ >12 are all underestimated. 
Better fits in the Case B regime are not obtained even if we 



' To derive the Cas e B pre dictions, the Fortran program provided 
by IStorev & Hummerl il995h and the required data files available 
at |http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=Vi764 have been 
used. 
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Fig. 4. Left panel: Comparison between the DK Cha Brackett decrement and a model assuming an expanding wind in LTE at 
T=3500K. The other parameters of the model are: r,=l R©, ro,,r=5 R©, ng"'=10'^cm"^, Vo=20kms"', V,„av=300kms"'and a=3 (see 
text for details). Right panel: Comparison between the DK Cha Paschen decrement and the same model assuming an expanding 
wind in LTE at T=3500 K and r,=3 R© (dashed-dot black line) and at T=2000 K and r,=2 R© (continuous red line). The other model 
parameters are the sa me as in the Brack ett decrement case. For comparison the weighted mean line ratios for the Paschen series on 
the CTTSs sample of iBarv et al.l (l2008l open green squares) were also included. 



Table 3. Brackett and Paschen line fluxes. 



Line 


FW • " 
(A) 


Fi 

(10"'^ erg s"' cm"^) 


PaS 


-12.9 


9.33 


Pay 


-18.6 


13.4 


Pa^ 


-23.0 


16.3 


Br 20 


-3.0 


2.04 


Br 19 


-1.5 


1.02 


Br 18 


-1.9 


1.30 


Br 17 


-2.0 


1.33 


Br 16 


-3.4 


2.29 


Br 15 


-3.4 


2.31 


Br 13 


-6.7 


4.56 


Br 11 


-8.0 


5.44 


Bry 


-11.2 


7.46 



" Equivalent width of the circumstellar H i line. Uncertainties are es- 
timated to be ± 1 A. 
* Fluxes corrected by an extinction of Ai/=11 mag 



increase the possible range of electron densities and tempera- 
tures: this suggests that the different behaviour of the Brackett 
lines in DK Cha with respect to the iBarv et aTl (12008^ sample 
of CTTSs is not caused by the higher luminosity/mass of the 
source. Unfortunately, we cannot prove which the behaviour of 
the Paschen lines at higher n„y, is, because only the Pa/3 Pay 
and Pa^ lines were detected at a sufliciently high S/N in a clear 
atmospheric region and not blended with other emission lines. 

Evidence that the H i lines in DK Cha cannot be reproduced 
by the optically thin regime assumed in the Case B recombi- 
nation is also provided by the observed (extinction corrected) 
PayS/Bry ratio. This ratio is 2.18+0.04, while for Case B it is 
never below 2.5: indeed, very low values of Pa/J/Bry are obtained 
only whe n the two lines are optically thick in LTE at low tem- 
peratures (iGatti et al.ll2006l;[Antoniucci et al.l201 lb . An estimate 
of the temperature and projected emitting area that gives rise to 



this thick emission can be obtained from the Pa/SfBry ratio and 
the Pa^ flux. Assuming a line-width of 150 km s"' , we derive 
T~3500K and a projected area of the order of a few solar radii. 
This indicates that the emission originates from a quite com- 
pact region with gas at low temperature. To verify the consis- 
tency of our previous statement, we also plot in Fig.|3](left panel) 
the expected Brackett decrement for a black-body at T=3500 K 
(red dotted dashed line), i.e., assuming that all Brackett lines are 
optically thick and at the same temperature. Evidently, the pre- 
dicted high-n„p lines are largely overestimated with respect to 
the Bry. Summarising this analysis, we conclude that the DK 
Cha Brackett decrement is consistent with emission at a quite 
low temperature where the lines are optically thick for the tran- 
sitions at low n,ip and optically thin for the higher n„p values. On 
the other hand, no additional information can be retrieved from 
our Paschen decrem ent (Fig.[3l r ight p anel). 

As discussed bv lBarv et al.l (l2008l) . the low temperature de- 
rived for the excitation of the Brackett and Paschen lines is not 
consistent with an origin in magneto-spheric accretion columns, 
where the temperatures are predi c ted to be in the r ange 6000- 
12000 K dMuzeroUe et al.lll998ah . iBarv et all (l2008h suggested 
that the direct absorption of high energy photons from the hot 
corona and/or accretion shock is potentially able to ionise the 
gas and produce the intense H i emission even at low tempera- 
tures. Qualitatively this seems applicable to the DK Cha case, 
which is a well known X-ray source with an X -ray luminosity of 
Lx -1.3x10^' ergs ' (iHamaguchi et al.ll2005h . 

4.2.2. Spherically symmetric expanding wind model 

Finally, to explore a different excitation scenario, we also com- 
pared the observed Brackett and Paschen decrements with the 
predictions of a simple model of an expanding wind, which was 
able to consistently r eproduce the H i li ne emission in the Class 
I source HHIOO-IRS (kisini et al.ll2004 . This model consists of 
a spherical envelope of ionised hydrogen in LTE , where the gas 
is expanding from an inner radius r,- to an outer radius Tout fol- 
lowing a velocity law of the type V(r)=Vo-i-(V„,a^--Vo)(l-(r,/r)"'), 
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where Vq is the initial velocity at the ba se of the wind and Vmax is 
the maximum velocity of the wind (see'N isini et alJl2004ll995l 
for more details). Under these conditions, the line decrement 
mainly depends on the optical depth of the line, which is influ- 
enced by the velocity law, the emitting region size, and the elec- 
tron density assumed at the wind base. Assuming a T~3500 K 
and an emitting region size with radius ~ 5 Rq, we are able 
to find quite a good agreement with all observed Brackett line 
ratios with a model having a Nmax ~ 300 kms"' and n"" ~ 
lO^^cm"^ (Fig.H] left panel). This model is also able to repro- 
duce the Pay6/Bry ratio as well as the extinction-corrected Pa/S 
flux of ~lxlO"" erg cm"^ s"'. 

We point out that although this is a toy model, the main re- 
sults do not depend much on the given assumptions. Although 
the main approximation here is spherical geometry, the high gas 
speed and the consequent high velocity gradients imply that the 
radiative interaction occurs only in the confined region of line 
formation. Consequently, the line emissivity depends only on lo- 
cal physical quantities and line ratios are not much affected by 
geometrical effects. 

The assumption of spherical wind is instead affecting the 
mass loss rate implied by our model: indeed, an estimate of M^. 
can be given from the continuity equation: 



Mw = AnrfVon" 



(1) 



Substituting the relevant parameters, we derive Mw - 
S.OxlO^^Moyr"'. Such a value is of the same order of mag- 
nitude as the DK Cha mass accretion rate derived in Sect. l4.4l 
However, Mw (eq.[T]) scales with the total surface area and could 
be a small fraction of the derived value if the emission comes, 
e.g., from a collimated wind. 

It is also worth noting that the same wind parameters used to 
fit the Brackett ratios do not reproduce the Paschen decrement, 
however (see, Fig.|4] right panel; black dashed line). Lower tem- 
peratures (T~2000K) or a greater envelope thickness (around a 
factor of 2-3) are required to match the previous wind solution. 
This may indicate that Paschen lines are emitted in a region that 
is more extended than the Brackett lines and/or that a temper- 
ature gradient is required. However, our model is probably too 
si mplistic to draw fu r ther c onclusions. 

iKwan & Fisched (1201 ih performed local excitation calcula- 
tions to obtain hydrogen line opacities and emissivity ratios, 
finding out that the Hi line emission mainly arises from a 
highly clumpy radial outflow. Following this work, the observed 
Pay/PayS ratio in DK Cha (Pa7/Paj6=0.82+0.07) is consistent 
with Pa/3 and Pay lines excited in a region with a total den- 
sity of roughly Nh ~L6x10"-6x10" cm"^ and an optical depth 
T/'fly=33.4-210, with both lines showing different excitation tem- 
peratures in the range 3240-3920 K and 3020-3940 K, respec- 
tively. This result resembles that from our simple model fairly 
well. 



4.3. Other detected lines. 

Previous observations of DK Cha have already shown some 
of the prope rt ies of its circumstellar material. For instance, 
iHughes et al.l (Il99lh reported the presence of a wind from 
this sou rce as shown by blue-sh ifted forbidden line emission, 
whereas [van Kempen et al.l (l2009l) detected a CO outflow. In our 
spectra, the presence of faint H2 emission together with the CO 
band-heads may also indicate the presence of a molecular out- 
flow. However, the near-IR H2 and CO emission could also be 



related to emission fro m the interaction of a wind with the sur- 
rounded medium (e.g., 'Garcia Lopez et al."2008') or/and emis- 
sion from the inn er warm circumstellar disk CCarr et al.,,1993t 
lNaiitaetalJ[l996l) . 



The optical and NIR spectrum shows several forbidden emis- 
sion lines that are typical of atomic jets, such as the sev- 
eral [Sii] [Oi] and [Ni] transitions between 0.6 and lyum. 
The pre sence of a jet fr om DK Cha was already pointed out 
by Hug hes et al.l ([l99Tb . who resolved [Sii] extended emis- 
sion showing significant velocity gradients. From the luminos- 
ity of the [Oi]0.630//m it is possible to estimate the mass 
flu x rate from the jet (M je,), following the relationship given 
bv ' HartiganetaD(ll995h . The main uncertainty using this ex- 
pression is caused by the assumed extinction: the Ay de- 
rived in Sect. l4.1l towards the central source can be consid- 
ered an upper limit to the extinction towards the region of 
the emission of the forbidden lines, which is likely to be lo- 
cated farther out in a less embedded environment. Therefore, 
we can only give an upper limit to the mass flux by assuming 
Av=ll mag for the extinction towards the forbidden line emis- 
sion region. In this way, we infer that Mje, is <10"^ yr"', 
assuming a distance to the source of ~178pc, an electron den- 
sity of -lO^cm"-' (because both emission from [Oi] and [Sii] 
have been detected, see Hartigan et al. 1995), and a veloc- 
ity in the plane of the sky ~32 kms~' (assuming i~ 18° 
and a radial veloci ty of — 98kms . Ivan Kempen etal1l20I0l: 
Hu ghes et al.lll99i]) . To better constrain the Mjet value in DK 
Cha, the [O i] 63 yum line has been used as well. This line has 
been detected in both ISO-LWS and HERSCHEL-PACS obser- 
vation s of DK Cha dvan Kempen et al.l I2OIOI: iLorenzetti et al.l 
1999!). Although this line is not alTected by extinction, the Mjei 
value derived from it represents an upper limit to the total 
mass transported along the jet because part of the [O i] 63 yum 
line could be excited in the PDR of the star (see e.g., 
ILorenzetti et al.l Il999l) . Assuming a [Oil 63 fim line flux of 
3.7xl0"'^ergs ' cm"^ (van Kempen et al. 1201 Ol) and the ex- 
pression reported in Cabrit (2002) to estimate Mjet, we found 
an upper hmit of ~3.6xlO"^ Mq yr"'. 



HERSCHEL-PACS observations of this source also revealed 
the presence of several H2O and OH lines interpreted as possi- 
bly generated from UV-heating of the material along the cavity 
walls and thus poin ting out the presence of a strong UV field 
from the protostar (Ivan Kempen et al] l2010l) . The optical and 
NIR spectra of DK Cha presented in this work show, indeed, sev- 
eral permitted transitions that are likely excited by UV pumping. 
The 1 0.845 yum and O i L 129 yum lines are two good examples: 
the presence of the bright Oi 1.129yum emission together with 
the absence of the Oi 1.132yum line indicate Hi Lyy6 fluores- 
cence of the 3d^D level of O i as the pumping mechanism for the 
1 0.845 yiim line. This suggests in turn that the H a line must be 
optically thick in this source because a large populatio n of the 
n=3 level of hydrogen is needed for this process to occur ( Grandil 
1975, 1980). Oxygen I fluorescence lines have been reported be- 
fore in_odierHerbig AeBe stars, as MWC349 and RMon (see 
e.g. iKellv et al.l[T994i) . In addition to the Oi lines, several Ni 
transitions at 1.013, 1.051 and 1.060/im, whi ch may be also 
excited by fluorescence, have been observed. IWalmslev et al.l 
(2000) have already reported the presence of some of these lines 
in the spectrum of the Orion bar They attributed the excitation 
of these lines to fluorescence within the ionisation front. 
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4.4. Accretion luminosity and mass accretion rate. 



Table 4. DKCha stellar parameters. 



The accretion luminosity (L„ce) in embedded YSOs is mostly 
derived from the luminosity of infrared lines such as the Bry 
and Pa/3 lines (Muzerolle et al. 1998b). The Bry line has been 
successfully used to derive the accretion luminosi ty in objects 
ranging from brown-dwarfs to s everal solar masses ( Calvet et al.l 
l2004t lG arcia Lopez et al.ir2006l) . however, the Pa /? line has been 
mostly used to compute the accretion luminosity in low-mass 
protostars (Gatti et al. 2008; Natta et al. 2006). The simultane- 
ous acquisition of the SOFl blue and red grism spectra allows 
us to derive Lace from the luminosity of the Pa/3 and Bry lines, 
which gives us the opportunity of comparing the inferred La^c 
values without introducing any uncertainty caused by intrinsic 
flux line variations of the source. 

With this aim, the empirical relations rel ating the accretion 
luminos ity with the Br y and Pa /3 luminosity dCalvet et al.ll2004| 
and Cal vet et al ] l2000h were used: 



log Lace I = 2.9 + 0.9 logLiBry) 
log Lace/ Lq = 2.80 + 1.03 logL(Pa/3). 



(2) 
(3) 



The luminosities of the Bry and Pa/3 lines (L(Bry), L(PayS)) 
were derived from the observed equivalent widths of the Br y and 
Pa/3 lines (EW(Bry), EW(Pay6)) once they were corrected from 
the intrinsic photospheric absorption contribution (see Sect l4.2l i. 

Using this procedure, we found a very good agreement be- 
tween the accretion luminosity computed from the Bry and Pa/3 
luminosities (see, Table|5j, with an average value of L„ef ~9 L©. 

Additional relationships, based on the luminosity of optical 
lines, have also been suggested for the determination of the ac - 
cretion luminosit y. For instan ce, Herczeg & Hillenbrand! (l2008h . 
[Oahm (2008) and lFang et al.l (l2009l) . derived relationships based 
on the luminosity of the [Oi], Can and Ho- lines, respectively. 
When we apply these relationships to the lines observed in 
DK Cha, we derive accretion luminosities much higher than the 
value estimated from the IR Hi lines, however (Lace ~104, 24 
and 107 Lq for the [Oi], Can and Ha lines, respectively), and, 
more important, much higher than the bolometric luminosity of 
this source. The origin of this discrepancy could be related to 
a wrong extinction estimate. For instance, decreasing the Ay of 
2mag, new Lgcc values of ~7.4, 4.5, 11.7, 18.6 and I2L0 are 
found for the Br y, Pay6, Ca 11, O i and i\a, respectively. Although 
these values are now lower than hboi, the accretion luminosities 
derived from the Bry and Pa/3 lines now disagree. Similar re- 
sults to the ones found in DK Cha (i.e. optical lines giving rise 
to Lflcc values higher than IR lines) have also been reported in 
a large sample of sou rces showing low extinc tion values located 
in the Chal/11 clouds (lAntoniucci et al.ll201 ll) . 

It is more likely, then, that the empirical relationships found 
for samples of low or very low mass objects cannot be extrap- 
olated to sources with larger mass and luminosity. In the high- 
luminosity objects different excitation mechanisms in addition 
to accretion can indeed be responsible for enhanced emission 
of the optical lines such as high chromosphere activity or direct 
excitation from stellar UV photons. 

From the accretion luminosity derived from the IR lines 
and stellar parameters shown in Table|4] a mass accretion rate 
(Mace) of ~3.5x lO~^Moyr~' was found from the expression 
dGullbring et al.l[T998l) 



ST 


D 

(pc) 


(K) 


(lI) 


Lfao/ 
(Lo) 


R* 
(Ro) 


(Mo) 


FO^ 


178 


7200^= 


18.62" 


29.4* 


2.77° 


2" 



" From:Spezzi et al. (2008) 

* From lvan Kempen et aU C2010h 



Table 5. Accretion properties. 



Ay 


L„cc(Pay6) Ueei^ry) 


Mace 


(mag) 


(Lo) 


(10-^Moyr-i) 


11.2 


9.5 9.1 


3.5 



Ma. 



LaccR* L _ R* 

GM, ~R~ 



(4) 



where R, and M» are the stellar radius and mass, G is the 
universal constant and a R„,=5R0 is the inner disk radius. 
The inferred Mace value is within the range of those found 
in classical Herbig Ae stars and intermediate mass T Tauri 
stars (IMTTS) where Mace in the range lO '^-lO'*^ Mp yr'' are 
typica lly found (D onehew & Brittain 201 ll iGarcia Lopez et all 
2006i: ICalvet et al.li2004 : Rodgers .i200lh . Interestingly, when one 
compares sources of roughly the same mass, the Mace value 
fou nd in DK Cha is hi gher th an that found in the IMTTS sample 
of Calvet et all (l2004 and in iGarcia Lopez etal.l (|2006') (where 
Mace values ~10~^ Mq yr"' were f ound) but lower than th a t com- 
puted for the Class I sample of IWhite & Hillenbrand! (|2004|) 
(Mace -10"'' Mo yr"'). This result is consistent with DK Cha be- 
ing in an evolutionary transition phase between a Class I and II 
source. 



5. Conclusions 

We have presented low-resolution spectroscopic observations 
of the embedded Herbig Ae star DK Cha covering the range 
0.6jL/m-2.4/vm. Its nearly face-on configuration has allowed us 
to have direct access to the star-disk system and to detect numer- 
ous emission lines tracing different emission regions and excita- 
tion conditions in the nearby surroundings of the protostar The 
extremely rich emission line spectrum makes this source very 
similar to low-mass very active T Tauri stars. We summarise our 
results as follows: 

- The DK Cha spectrum shows several permitted (e.g., 
Ca ii,Mg i,Na i), forbidden (e.g., [Fen], [Sii]) and molecular 
emission lines (CO and H2). The Hi lines are those more 
numerous along the spectrum. In addition, some of the per- 
mitted lines, such as the 1 0.845 fim line, were identified as 
being excited by fluorescence. 

- To constrain the origin of the H i lines, Brackett decrement 
plots were constructed and compared with different excita- 
tion mechanisms. Case B recombination is not able to fit all 
observed Brackett ratios, while optically thick emission in 
LTE at a temperature of ~ 3500 K is required to account 
for the observed PayS/Bry ratio. A simple model for an ex- 
panding gas is able to reproduce all observed Brackett ra- 
tios by assuming that this gas is located very close to the 
central source and has a high density at its base, of the or- 
der of lO'^cm"^. The same parameters used to reproduce the 
Brackett ratios cannot account for the Paschen line ratios. 

- From the Pa/3 and Bry line luminosity we de rived the ac- 
cretion luminosity using the expressions from ICalvet et~an 
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(1200(1 12004 . We found a very good agreement between 
the Lace derived from both expressions, with an aver- 
age value of La„~9L0. From the accretion luminosity 
a Mace~3.5xl0"^Moyr"' was found. The measured Mace 
value is consistent with DK Cha being in an evolutionary 
transition phase between a Class I and II source. 
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